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面向 B 超 图 像 分 割 的 动态 权重 因子 水 平 集 方法 
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摘要 :目的 研究 和 改进 水 平 集 方法 ,实现 B 超 图 像 中 病灶 区 域 的 准确 快速 分 割 。 方 法 分 析 已 有 水 平 集 方法 在 B 超 图 像 处 理 中 
的 局 限 ,基于 区 域 水 平 集 的 优点 ,将 信息 论 中 的 炉 引 入 图 像 处 理 ,定义 动态 权重 因子 ,准确 反映 局 部 灰 度 阶梯 变化 状况 ,定量 度 
量 轮廓 线 像素 点 分 别 受到 趋向 目标 .背景 区 域 的 两 种 作用 力 的 动态 权重 ,将 其 融合 到 区 域 水 平 集中 ,迭代 引导 曲线 形变 和 位 
移 。 由 于 B 超 图 像 病灶 分 割 属于 指定 区 域 的 局 部 分 割 , 所 以 将 计算 约束 到 局 部 范围 ,从 而 明显 降低 运算 代价 。 结 果 动态 权重 


因子 水 平 集 方法 能 够 较 好 分 割 B 超 图 像 中 的 病灶 区 域 ,与 几 种 主流 水 平 集 方法 相 比 ,本 文 方法 精度 更 高 ,时 间 复 杂 度 更 小 。 结 
论 动态 权重 因子 方法 能 够 更 合理 准确 地 判断 病灶 边界 像素 点 ,局 部 计算 策略 有 效 地 提高 了 分 割 效 率 。 
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Abstract: Objective To modify the level set method for precise and fast segmentation of B-type ultrasound image lesions. 
Methods Based on the best of region level set method, entropy in the information theory was introduced into image 


processing to define a dynamic weighting factor that responded to the gradient change of the local gray levels to evaluate the 
dynamic degree of driven force on each pixel on the contour to the target and background areas. The dynamic weighting 
factors were reconciled into the regional level set model and led the contour to deform and move during the iterations. As 


lesion segmentation was classified as local segmentation of a specific area, the calculation was restrained to the local sphere 
abide by the contour, which greatly reduced the calculation complex. Results Experiments on B-type ultrasound images 
showed improved results of the proposed method with a better accuracy and less time consumption compared with several 
current level set methods. Conclusion Level set method reconciled with dynamic weighting factor allows a better evaluation of 


the lesion contour pixels, and the local calculation strategy results in an enhanced segmentation efficiency. 
Key words: level set method; dynamic weighting factor; image entropy; B-mode ultrasound image extraction 


B 型 超声 对 液体 与 实质 性 组 织 表 现 出 显著 的 图 像 
差别 ,在 发 现 肿 物 .畸形 结石 等 引致 局 部 结构 明显 形态 
改变 的 疾病 方面 有 着 明显 优势 ,无 创 方便 廉价 实时 等 优 
点 使 得 它 在 体检 和 医疗 诊治 中 的 应 用 十 分 广泛 "。 由 
于 成 像 原理 和 过 程 中 的 多 种 因素 影响 ,B 超 图 像 常 表现 
出 对 比 度 较 低 .噪声 干扰 严重 边界 模糊 等 缺点 ,自动 分 
割 在 B 超 中 的 应 用 具有 较 大 局 限 性 。 人 们 从 多 种 角度 
采用 多 种 方法 进行 B 超 图 像 的 自动 分 割 探索 ,如 区 域 增 
长 法 .动态 国 值 分 割 法 模糊 聚 类 法 、 分 水 岭 法 .HMM 融 
合法 EM 算法 和 高 斯 混合 分 布 模型 ,小波 变换 法 等 , 现 
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有 解决 方案 尚 存在 较 大 的 改进 空间 。 基 于 变 分 的 水 平 
集 方法 是 近年 来 发 展 迅 速 的 一 个 分 支 , 它 将 图 像 分 割 问 
题 表 达 为 能 量 函 数 的 最 小 化 问题 ,由 变 分 原理 将 其 转化 
为 偏 微分 方程 的 求解 。 相 比 于 传统 的 区 域 分 割 方法 ， 
水 平 集 方法 允许 在 能 量 函 数 中 加 入 与 图 像 内 容 有 关 
的 约束 条 件 等 综合 信息 ,获得 更 加 准确 和 合理 的 分 
HAR 

水 平 集 方法 的 主要 思想 是 将 移动 的 界面 作为 零 水 
平 集 胖 入 高 一 维 的 水 平 集 函 数 中 ,由 封闭 超 曲 面 零 水 平 
集 的 演化 方程 得 到 水 平 集 函 数 的 演化 方程 ,通过 水 平 集 
函数 曲面 的 演化 来 隐 含 地 求解 轮 廊 曲线。 水平 集 方法 
可 分 为 全 局 水 平 集 ;边界 水 平 集 和 区 域 水 平 集 三 类 。 

以 CV 模 型 "为 代表 的 全 局 水 平 集 模型 利用 图 像 的 
全 局 灰 度 信息 驱动 曲线 形变 和 移动 ,主要 优点 是 当 两 种 
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区 域 的 灰 度 差别 较 大 时 具有 较 强 的 抗 噪 能 力 ,能 够 有 效 
分 割 出 模糊 边缘 ;对 初始 轮廓 线 位 置 不 敏感 ;经 理论 证 
明 能 够 收敛 到 全 局 最 优 。 但 是 , 它 以 图 像 灰 度 分 布 均 
^] ,前景 背 景 有 明显 差异 为 前 提 条 件 , 不 能 直接 用 于 了 B 
超 图 像 的 病灶 分 割 。 以 GAC 模 型 2 为 代表 的 边界 模型 
利用 目标 边界 的 梯度 变化 信息 引导 曲线 运行 ,优点 是 曲 
线 在 平滑 区 域 演化 很 快 ,而 在 目标 边界 处 停止 运动 ,在 
低 噪 图 像 中 分 割 准 确 且 运算 代价 小 ,但 在 高 噪 情况 下 容 
易 受 到 噪声 干扰 而 使 曲线 停顿 在 错误 的 位 置 上 。 以 
RSF 模 型 "为 代表 的 区 域 水 平 集 模 型 利用 曲线 内 外 部 局 
部 区 域 的 拟 合 能 量 信息 促使 曲线 形变 ,能 够 较 好 地 处 理 
灰 度 不 匀 的 情况 ,但 运算 代价 比 全 局 和 局 部 模型 显著 增 
加 ,使 用 的 核 函 数 的 方差 为 定 值 ,容易 陷入 局 部 极 值 
解 。 近 年 来 围绕 RSF 模 型 的 改进 主要 分 为 两 个 方向 : 核 
函数 的 选择 ,如 用 其 他 函数 (均值 方差) 来 代替 高 斯 分 
布 函 数 呈 ;融合 其 他 方法 ,如 文献 [6-7] 向 RSF 模 型 中 加 
入 了 模糊 聚 类 方法 ,在 一 定 程度 上 提高 了 像素 点 的 判断 
精度 ,文献 [7] 还 提出 了 快速 计算 方案 ,但 对 复杂 图 像 仍 
存在 模糊 度 选取 问题 ,可 能 造成 过 分 割 或 从 分 割 。 
文献 [8-10] 融 合 了 CV 和 RSF 模 型 ,可 利用 两 者 的 优 
点 ,但 双 模型 的 计算 量 显 著 增加 。 

本 文 针 对 B 超 图 像 特点 ,利用 区 域 水 平 集 模型 的 优 
点 ,定义 反映 曲线 像素 点 受到 目标 和 背景 两 个 区 域 动态 
权重 的 因子 ,通过 概率 分 布 估计 模型 计算 像素 点 被 动态 
权重 的 程度 ,以 此 为 依据 对 像素 点 进行 区 域 划分 ,以 提 
高 水 平 集 曲 线 抗 干扰 能 力 和 定位 模糊 边界 能 力 ; 同 时 本 
文 将 B 超 图 像 目 标 分 割 看 作对 感 兴趣 目标 区 域 的 局 部 
分 割 ,将 水 平 集 的 计算 求解 约束 到 局 部 分 割 范围 ,减少 
计算 量 ,实现 特定 目标 区 域 的 高 效 提 取 。 


1 动态 权重 因子 的 定义 和 融入 区 域 水 平 集 模 型 

设 C={ C, . C, 是 岁 像 [关于 轮廓 线 C 的 两 个 划 
分 ,C, 代表 轮廓 线 内 部 即 目标 区 域 ,C， 代 表 轮 廓 线 外 
部 即 背 景区 域 ,RSF 模 型 的 能 量 函 数 定 义 为 : 

Est Ao e)" > h &c-»| I()-fo[ ay. i-1,2 (1) 


其 中 ,9 和 Q, 分 别 为 轮 廊 线 内 部 和 外 部 区 域 ,K 为 
核 函 数 ,一 般 选 择 为 高 斯 函数 ,f(x) 和 户 (X?) 为 像素 点 x 处 
的 拟 合 值 ,4 和 ,为 常 系数 。 

分 析 图 像 ,在 目标 内 部 和 背景 内 部 ,局 部 灰 度 分 布 
一 般 是 均匀 的 ,而 在 目标 与 背景 的 交界 处 ,局 部 灰 度 分 
布 的 差异 性 较为 显著 。 因 此 , 某 个 区 域内 部 的 灰 度 分 布 
越 均 衡 , 则 越 应 当 作 为 一 个 整体 被 划分 , 即 对 曲线 像素 
点 的 权重 作用 越 大 ,曲线 最 后 应 当 包 绕 在 该 区 域 的 边界 
上 ;反之 ,如 果 某 个 区 域内 部 的 灰 度 分 布 不 均衡 即 差异 
较 大 , 则 表明 它 很 可 能 既 包 含 了 目标 部 分 又 包含 了 背景 


部 分 ,或 者 带 有 强 噪声 点 ,曲线 不 应 当 收 拢 到 该 区 域 的 
边界 上 ,也 就 是 该 区 域 对 像素 点 的 权重 作用 应 当 较 小 。 
RSEF 模 型 中 轮廓 线 内 部 和 外 部 区 域 的 作用 力 的 权重 由 
常 系数 0 和 41 来 规定 ,一 则 需要 手工 指定 ,依赖 经 验 值 ， 
二 则 不 利于 动态 调整 。 图 像 区 域 灰 度 分 布 均衡 度 常用 
的 度量 指标 有 均值 方差 . 炉 等 ,其 中 烂 (entropy) 作 为 概 
率 系统 的 不 确定 程度 指标 , 当 系 统 越 确定 时 , 焙 值 越 小 ， 
反之 ,系统 越 不 确定 , 炉 值 越 大 。 在 统计 学 中 , 炉 的 玻 尔 


S-x-LnQ (2) 
其 中 ,x 为 ( 玻 尔 兹 曼 ) 常 数 ,Ln 为 对 数 运 算 符 ,Q 为 系统 分 
子 的 状态 数 。 


炉 的 概念 同样 可 以 映射 到 图 像 灰 度 分 布 的 均衡 性 
上 。 图 像 中 一 个 区 域 的 灰 度 分 布 越 均匀 , DUE EGER], 
Bu RUREK. Ha iie fe ee rp fee RR RE RE, EE 
Beo LACE SIMA AEBEC PASA — BEA E , SLATE RIE e UR 
图 像 中 的 几何 形状 .像素 的 空间 分 布 等 重要 的 处 理 
信息 ”””。 本 文 利用 业 的 基本 含义 ,对 它 在 图 像 处 理 中 
的 定义 做 了 修改 :图 像 的 箭 表 示 为 图 像 灰 度 级 集合 的 比 
特 平 均 数 ,描述 了 图 像 信 源 的 平均 信息 量 , 即 图 像 的 信 
息 丰 富 程 度 。 设 x(i,j) 为 图 像 的 像素 点 ,图 像 炉 对 应 着 
图 像 分 布 的 直方 图 ,近似 于 采样 均值 。 全 局 采样 均值 不 
能 准确 反映 图 像 各 个 部 分 的 像素 分 布 情况 和 变化 情况 ， 
所 以 ,本 文 设计 了 局 部 烂 , 即 局 部 区 域 采 样 均值 ,局 部 业 
的 定义 为 


H(p-- > pj hnp,G.p(2- 
ze B(x, e) 


|(z € B(x, e:1(x)=z}| (3) 
LOI 


其 中 ,p,(z) 束 示 以 像素 点 x 为 中 心 的 半径 为 的 窗口 中 的 
灰 度 级 概率 密度 。 局 部 炉 比 均值 和 方差 能 更 有 效 地 反 
映 灰 度 阶梯 变化 状况 。 轮 廓 线 上 的 像素 点 x 受到 局 部 焙 
的 作用 ,用 一 个 新 的 因子 :动态 权重 因子 来 表示 。 动 态 
权重 因子 定义 为 半径 为 e 的 区 域 的 局 部 炉 A.(7(0)): 
AATO)=— 2 pli, DInpli,)) (4) 


其 中 Q 为 关于 像素 点 x 的 半径 为 e 的 区 域 , /为 像素 点 zx 
的 坐标 。 
RSF 模 型 的 水 平 集 表示 为 
2 
Eg (9 f 09-509) 7 2- [, KG) -AO MpONdy, i=1,2 


(5) 
其 中 Mi(o) 7 H(9).M;() - 1— H (9). Hx Heaviside PR 
SGH o)- Hi = arctan(*) 
fG9- KOM POI G)] ] 
i K(x) M (9(x)) 
KGAA RSF EC BE E KA, HH T RRR 
TARRE, FEAE NT DK PRERA: 
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2 
EC. f 09. f,09) E= 4C D! [. G, 634,009) I0) -/,c0 | dy 
k=1 k 
- Yo GD Y, pi. Dnpi, DO-A dy 
k=1 s ijeQ, 


(6) 
其 中 ,G.(7(x)) 为 以 像素 点 x 为 中 心 半径 为 e 的 区 域内 的 高 
斯 分 布 ,Q 和 @ ,分别 为 轮廓 线 内 部 和 外 部 的 区 域 。 

可 以 看 到 ,(6) 式 中 的 两 项 分 别 为 像素 点 受到 曲线 
内 部 和 外 部 区 域 的 动态 权重 力 , 方 向 相反 , 当 二 力 相 等 
时 (6) 式 能 量 汤 数 E, 达 到 最 小 值 ,此 时 水 平 集 停止 在 目 
标 和 背景 的 边界 上 ,得 到 分 割 轮廓 线 。 动 态 权 重 因子 的 
作用 相当 于 一 个 变化 的 系数 , 随 着 轮廓 线形 变 和 位 移 到 
图 像 的 不 同 区 域 , 动 态 调节 两 项 力 的 大 小 ,促使 轮廓 线 
较 快 地 向 目标 边界 移动 和 形变 。 

引入 水 平 集 表示 法 ,得 到 动态 权重 因子 水 平 集 演化 
方程 : 

-ó(g) Soy R | p. jn p. DIO -AW dy 


k=1 


(7) 
其 中 5,(p)=H,0)= 上 22s。 这 样 就 把 图 像 目标 分 守 
E X 


问题 转化 为 极 小 化 能 量 孙 数 (7) 的 计算 。 


2 计算 方法 

区 域 水 平 集 模型 将 图 像 中 的 所 有 目标 分 割 出 来 , 计 
算 量 很 大 。 本 文 根 据 B 超 图 像 病灶 的 分 割 要 求 , 将 初始 
轮 廊 设 为 一 个 接近 并 包 绕 B 超 图 像 目标 的 闭合 曲线 ,为 
简便 起 见 ,统一 采用 手工 勾勒 包围 目标 规则 图 形 (和 矩 形 、 
椭圆 形 或 圆 形 )。 为 提高 计算 效率 ,本 文采 用 局 部 计算 
方案 ,定义 基于 局 部 近似 符号 距离 函数 ,放宽 对 初始 水 
平 集 函数 的 要 求 , 不 再 要 求 初始 水 平 集 必须 是 符号 距离 
函数 。 如 果 直 接 对 轮廓 线 上 的 每 个 像素 点 都 做 对 数 及 
求 和 运算 ,求解 水 平 集 式 需要 离散 化 迭代 实现 求解 偏 微 
分 方程 ,计算 量 很 大 ,因此 采用 局 部 计算 思路 以 减少 运算 
量 和 运算 时 间 。 本 文采 用 改进 的 FTC(Fast two-cycle) 
算法 ,在 不 用 求解 偏 微分 方程 的 基础 上 近似 实现 
水 平 集 函 数 的 进化 ,定义 局 部 近似 符号 距离 也 数 
(LASDF) 为 


d, 点 z 为 外 部 点 

1, 点 z 为 轮廓 线 外 邻 域 点 
0， 点 x 在 轮廓 线 上 (8) 

-1,x 为 轮廓 线 内 邻 域 点 

-d , 点 xz 为 内 部 点 

LASDF 确 定 当前 像素 点 到 水 平 集 曲 线 的 最 短 距离 
和 该 距离 的 符号 (曲线 内 为 负 ,曲线 外 为 正 ) ,计算 过 程 
只 关注 包 绕 水 平 集 曲 线 的 一 条 半径 为 &( 宽 度 为 2d) 的 狭 
窄带 区 ,由 于 d<<n ,所 以 可 将 计算 复杂 度 由 00- In SWÈ 


qx, y- 


少 为 O(nd1n 5),n 为 图 像 中 的 像素 点 数 ,5 为 计算 每 一 个 


像素 点 到 水 平 集 曲 线 的 最 短 距 离 的 耗费 。 水 平 集 迭 代 
过 程 伪 代 码 如 下 : 


function phiy=activecontourCV(u0， center, d, 
side, d it, m it, m name) 

% 输入 图 像 为 double 型 , 灰 度 为 0 一 255 

% 水 平 集 初始 化 为 椭圆 形 ,center 为 轮廓 线 原 点 ,d 
为 窄带 半径 ,side =1 表 示 目 标 在 起 始 轮廓 线 内 ,=0 表 示 
在 外 

% 初始 化 参数 

ITERATIONS-2000; % 人 迭代 次 数 

delta_t=0.1; % 时 间 步 长 

% 轮 廓 内 外 能 量 参数 

lambdal=1; lambda2=1; nu=0; d=3; 

% 曲 率 项 参数 

h=1; h_sq=h^2; epsilonz1; mu=0.01 * 255^2; 

% 初始 化 LASDF 

phizinitLASDF(size(uO), center, d, side); 

for iiz1: ITERATIONS; 

多 显示 并 保存 图 像 

if(mod(ii-1, d it)2-0) 

segim-createim(uO, phi); 


imshow(segim); 

segim-createim(uO, phi); 

filename-strcat(m, name, sprintf( '%d', i, '.jpg' ); 

imwrite(segim, filename); 

end; 

dirac delta t-delta t * diracfunction(phi, epsilon); 

% 计算 轮廓 线 内 外 能 量 

[inside, outside]=calcenergyf(u0, phi, epsilon); 

energy term--nu -lambdal .* inside + lambda2.* 
outside; 

dx centralz(circshift(phi, [0, -1])-circshift(phi, [0, 
1py2; 

dy. centralz(circshift(phi, [-1, 0])-circshift(phi, [1, 
0]))/2; 

abs grad phi-(sqrt(dx central.^2 + dy_central.^2) 
4-0.00001); 

x-dx central ./ abs grad phi; 

yzdy. central ./ abs grad phi; 

grad term-(mu / h sq) .* divergence(x, y); 

phi-phi + dirac delta t .* — (grad term. + 
energy term ); 

phiy(:,:,11)=phi; 

end; 


% 检 测 零 水 平 集 的 位 置 
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function front=isfro(phi) 

[n, m]=size(phi); 

front=zeros(size( phi)); 

for i=2: n-1; 

for j22: m-1; 

max Valzmax(max(phi(i:i71, j:j-1))); 
minValzmin(min(phi( i:i+1, j:j-1))); 

front(i, j)2((maxVal*0) & (minVal «0)) | phi(i, j)- 


end 

end 

% 初 始 化 水 平 集 函数 ,局 部 近似 符号 距离 函数 
function phi=initLASDF (imsize, center, d, side) 
m-imsize(1); n=imsize(2); 
phi=zeros(imsize); 

for i=1:m; 

for j=1:n; 

distance=sqrt(sum((center-[ 1, j ] ).^2)); 
phiGi, j)2distance-d; 

if(side == 0) 

phi i, j )=-phi( i, j ); 

end 

end 

end 

% dirac 

function yzdiracfunction(x, epsilon) 
y=(1./ pi) .* (epsilon./(epsilon.^24x.^2)); 
% 求 轮廓 内 外 能 量 

function value = heavisidef(z, epsilon) 
value-0.5 .* (1+(2./pi) .* atan(z/epsilon)); 


function[inside, outside] -calcenergyf(uO, phi, 


epsilon) 
H  phi-heavisidef(phi, epsilon); 
H. phi. minusz1-heavisidef(phi, epsilon); 
A-calcuProb(d); 
cl=A.*sum(sum(u0.* H. phi))/ sum(sum(H . phi)); 
Cc2-A.*sum(sum(u0.*H, phi minus))/sum(sum 
(H. phi, minus)); 
insidez(u0-c1).^2; 
outside-(u0-c2).^2; 


3 结果 

分 割 实验 待 测 图 像 来 自 南 方 医科 大 学 附属 医院 影 
像 库 , 通 过 对 相同 B 超 图 像样 本 进行 分 制 来 检验 本 文 模 
型 的 效果 ,并 与 几 种 水 平 集 模型 进行 比较 。 编 程 平台 为 
Matlab7, 计 算 机 配置 为 ntel@ Core™(2) Duo CPU. 3.0 
GHz,2 GB RAM, Windows XP 操作 系统 。 所 有 实验 采 
用 相同 的 终止 阔 值 人 已 = 10 ”或 最 大 迭代 次 数 2000( 考 
虑 到 用 户 等 待 时 间 )。 图 1~3 为 实验 结果 举例 ,比较 
CV 模型 RSFEUU SLGS EU" LGF FRANU Z 
文 模型 和 专业 医师 手工 分 割 结 果 。 所 有 模型 均 按 原 
文献 方法 进行 相应 参数 设置 ,其 中 公共 参数 :时 间 步 
长 At=0.1 ,空间 步 长 Ah=1, 长 度 惩罚 项 参数 
115 0.01 * 255^2~0.1 * 255^2, 邻 域 半径 se 根 据 图 像 内 容 
调整 为 最 合适 值 (为 5 或 7), 尽 可 能 达到 每 种 模型 最 优 
的 分 割 效果 。 本 文 方法 中 的 罕 带 半径 dg 取 3, 即 计算 水 平 
集 曲 线 内 外 各 3 个 像素 的 LASDE。 

图 1A 子 图 为 待 测 样本 390*290 像 素 的 前 列 腺 增生 
B 超 图 像 ,分 割 目 标 为 中 下 部 团 块 暗 区 。 图 2A 子 图 为 
待 测 样本 488*415 像 素 的 盆腔 (卵巢 ) 也 超 图 像 ,分 割 目 
标 为 中 下 部 的 卵黄 吉 部 分 。 图 3A 子 图 为 待 测 样本 
615*369 像 素 的 盆腔 (子宫 肌 瘤 妆 超 图 像 ,分 割 目标 为 


1 分 割 测 试 1 


Fig.1 Segmentation test 1. A: Sample image to be segmented; B: Initial contour by manual; C: Segmentation 


result by CV model; D: Segmentation result by RSF model; E: Segmentation result by SLGS model; F: 
Segmentation result by LGF model; G: Segmentation result by the proposed model; H: Result of manual 


segmentation. 
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图 2 分 割 测试 2 


Fig.2 Segmentation test 2. A: Sample image to be segmented; B: initial contour by manual; C: Segmentation 
result by CV model; D: Segmentation result by RSF model; E: Segmentation result by SLGS model; F: 
Segmentation result by LGF model; G: Segmentation result by the proposed model; H: Result of manual 


segmentation. 


图 3 分 割 测试 3 
Fig.3 Segmentation test 3. A: Sample image to be segmented; B: initial contour by manual; C: Segmentation 
result by CV model; D: Segmentation result by RSF model; E: Segmentation result by SLGS model; F: 
Segmentation result by LGF model; G: Segmentation result by the proposed model; H: Result of manual 


segmentation. 


中 部 偏 右 的 子宫 肌 瘤 体 。 图 4A 子 图 为 待 测 样本 361* 
252 像 素 的 胆 圳 B 超 图 像 ,分 割 目标 为 中 部 横向 暗色 团 
块 区 域 。 待 分 割 图 像 均 存 在 明显 的 灰 度 分 布 不 匀 ,存在 
不 同 数量 的 多 处 伪 影 ,存在 散布 的 不 同 强度 的 噪声 杂 
质 ,病灶 部 分 边界 比较 模糊 。 

根据 医学 诊断 特异 性 的 标准 采用 4 个 精度 指示 
器 : 真 阳 性 率 (TP) , 假 阳 性 率 (FP) , 假 阴 性 率 (FEN) , 真 阴 
性 率 (TN) ,结合 临床 诊断 要 求 定 义 4 个 分 割 精度 指标 : 

可 接受 率 accept rate: 由 专业 医师 判断 自动 分 割 结 
果 是 否 可 接受 ( 含 经 过 后 续 处 理 ) ,是 否 有 参考 价值 。 


测 准 度 Accuracy= TP n = TAI FN’ 反映 检验 
阳性 和 阴性 与 按 该 盘 检 实验 的 标准 值 的 比值 。 
灵敏 度 Sensirivipy = 7 人 ,又 称 真 阳性 率 ,反映 


入 检 实 验 发 现 病灶 的 能 力 。 


特异 度 Specificiy= z pp ,又 称 真 阴性 率 ,反映 
入 检 实 验 确定 非 病 灶 的 能 

表 1 为 5 种 模型 处 理 图 1~4 的 迭代 次 数 . 耗 时 和 分 
割 精度 指标 。 对 于 分 割 结果 存在 多 个 连续 区 域 的 , 取 最 
大 的 区 域 或 医师 手工 指定 的 区 域 。 分 制 失 败 案例 的 指 
标 无 法 计算 ,用 “" 表 示 。 

从 南方 医院 临床 医学 影像 图 像 库 中 随机 抽取 20 
幅 B 超 图 像 ,用 5 种 模型 分 别 进 行 处 理 ,以 获得 处 理 
精度 的 统计 结果 (只 取 全 部 分 割 成 功 的 结果 )。 由 于 
各 幅 图 像 大 小 .内 容 和 质量 不 一 致 ,所 以 不 进行 迁 代 
次 数 和 耗 时 的 比较 ,只 进行 处 理 结 果 精 度 的 统计 比 
较 ,结果 见 表 2。 


4 讨论 


从 以 上 实验 及 统计 数据 可 以 得 出 结论 :在 分 割 精 度 
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图 4 分 割 测试 4 

Fig.4 Segmentation test 4. A: Sample image to be segmented; B: initial contour by manual; C: Segmentation 
result by CV model; D: Segmentation result by RSF model; E: Segmentation result by SLGS model; F: 
Segmentation result by LGF model; G: Segmentation result by the proposed model; H: Result of manual 
segmentation. 


表 1 5 种 模型 的 迭代 次 数 、 耗 时 ( 秒 ) 和 分 割 精度 比较 


Tab.1 Comparison of 5 models in iteration times, time consumption (seconds) and segmentation precision 


Model Iteration times Time consumption Accuracy Sensitivity Specificity 
Test 1 
CV 2000 132.785 0.892 0.865 0.918 
RSF 640 54.207 0.921 0.942 0.901 
SLGS 80 42.431 0.915 0.929 0.901 
LGF 50 21.854 0.718 0.944 0.547 
The proposed model 230 32.866 0.944 0.949 0.925 
Test 2 
CV 390 58.695 0.915 0.898 0.933 
RSF 80 18.958 0.911 0.913 0.910 
SLGS 40 61.624 0.915 0.929 0.901 
LGF 15 19.890 0.779 0.879 0.715 
The proposed model 90 14.594 0.927 0.948 0.907 
Test 3 
CV 160 25.736 0.885 0.928 0.847 
RSF 90 23.547 0.917 0.915 0.919 
SLGS 40 67.783 - - - 
LGF 15 23.359 - - - 
The proposed model 100 17.612 0.940085 0.950 0.931 
Test 4 
CV 1900 107.053 0.860 0.916 0.804 
RSF 420 37.429 0.912 0.928 0.916 
SLGS 50 34.797 0.884 0.932 0.842 
LGF 30 17.375 0.805 0.879 0.747 
The proposed model 150 23.549 0.933 0.945 0.934 


表 2 5 种 模型 对 20 幅 随机 图 像 的 平均 分 割 精度 比较 


Tab.2 Comparison of average segmentation precision of 20 random images by 5 models 


Model Accept rate Average accuracy Average sensitivity Average specificity 
CV 100% 0.888 0.895 0.875 
RSF 100% 0.914 0.904 0.902 
SLGS 70% 0.878 0.902 0.890 
LGF 75% 0.741 0.891 0.670 


The proposed model 100% 0.937 0.949 0.926 
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方面 ,本 文 模型 不 同 程度 地 优 于 其 它 模型 ,各 次 实验 的 
可 接受 率 与 CV .PSF 模 型 相当 ,基本 上 可 全 部 被 人 工 接 
受 , 而 平均 测 准 度 .平均 灵敏 度 和 平均 特异 度 几 乎 都 高 
于 其 他 模型 ,其 他 模型 的 总 体 分 割 精 度 排 序 依次 是 : 
RSF 模 型 .CV 模 型.SLGS 模型 和 LGF 模 型 ,后 两 种 模 
型 效果 受 B 超 图 像 内 容 或 质量 的 影响 较 大 , 均 有 分 割 不 
成 功 的 样本 。 表 明 本 文 模型 较为 适合 分 割 B 超 病变 组 
织 , 同 时 得 检 非 病灶 的 能 力 较 强 。 分 析 原 因 ,CV 模 型 借 
助 全 局 信息 提高 了 抗 噪 能 力 ,但 未 考虑 局 部 区 域 信 息 ， 
水 平 集 曲 线 粗略 逼近 目标 边界 ,无 法 准确 分 割 灰 度 不 
^J 、 伪 影 干 扰 和 目标 背景 差别 小 的 目标 , 且 需 要 对 所 有 
图 像 点 进行 迭代 计算 , 耗 时 较 多 ;RSF 模 型 运用 局 部 拟 
合 信息 能 够 较 好 反映 区 域 图 像 属性 ,在 初始 轮廓 接近 目 
标 时 分 割 较为 准确 ,但 依靠 常数 权重 确定 形变 作用 的 比 
例 关系 ,不 能 准确 反映 曲线 位 于 不 同 区 域 时 的 动态 变 
化 ,对 模糊 边界 和 噪声 点 容易 出 现 错 误 ;SLGS 模 型 能 够 
较 好 地 提取 弱 边 界 ,又 有 一 定 的 抗 噪 性 ,但 由 于 采用 了 
全 局 分 割 策略 ,所 以 无 法 剔除 背景 中 与 病灶 相似 的 成 
分 ,导致 较 多 的 碎片 结果 和 误 判 断 ;LGE 模 型 通过 小 邻 
域 高 斯 函数 能 够 较 准 确 地 表现 局 部 的 不 规则 变化 ,但 容 
易 产 生 大 量 元 余 轮 廓 碎片 而 无 法 连 成 完整 区 域 ;本 文 方 
法 在 RSF 模 型 的 基础 上 引入 图 像 箭 的 作用 ,计算 动态 权 
重 因子 ,能 够 根据 曲线 不 同时 刻 的 变化 位 置 实时 调节 作 
用 力 的 权重 ,更 为 合理 地 判断 弱 边 界 上 像素 点 的 归属 ， 
并 且 具 有 一 定 的 抗 噪 能 力 ,分 割 结果 最 接近 专业 医师 的 
手工 分 割 结果 ,效果 最 好 。 

在 运算 代价 与 分 割 速度 方面 ,本 文 模型 迭代 次 数 多 
于 SLGS 模 型 和 LGF 模 型 , 少 于 CV 模型 和 RSF 模 型 ， 
实验 耗 时 在 部 分 情况 下 仅 比 LGF 模 型 略 多 , 比 其 他 模 
型 均 有 明显 减少 ,说 明 本 文 模型 的 计算 法 则 以 及 局 部 计 
算 方 案 具 有 一 定 的 降低 计算 复杂 度 的 能 力 ,有 效 地 减少 
了 总 计算 量 。 


5 结论 

本 文 提出 一 种 面向 B 超 图 像 病灶 分 割 的 动态 权重 
因子 区 域 水 平 集 方法 ,利用 图 像 炉 计算 动态 权重 因子 ， 
动态 更 新 轮廓 线 内 外 区 域 对 像素 点 的 动态 权重 作用 , D 
定 曲线 移动 方向 和 形变 趋向 。 在 计算 方法 上 ,不 采用 传 
统 的 解 偏 微分 方程 方法 ,不 需要 扫描 整个 图 像 , 仅 对 轮 
廓 曲线 内 外 邻 域 点 实行 转换 判断 处 理 , 更 新 水 平 集 演化 
趋势 , 减 小 了 计算 量 。 与 几 种 水 平 集 模型 的 比较 ,本 文 
模型 在 平均 测 准 度 平均 灵敏 度 和 平均 特异 度 方面 都 表 
现 较 好 ,并 且 起 伏 很 小 ,不 容易 受 图 像 内 容 和 质量 的 影 
啊 , 较 好 地 解决 了 B 超 图 像 灰 度 分 布 不 勺 、. 低 对 比 度 和 


模糊 边界 的 问题 。 本 文 方法 的 局 限 性 是 动态 权重 因子 
是 依据 区 域 灰 度 分 布 进行 计算 的 ,而 人 眼 视 觉 注 意 机 制 
的 研究 表明 ,对 人 了 眼影 响 较 大 的 低层 视觉 显著 性 特征 还 
包括 亮度 ,对比度 .目标 区 域 大 小 .目标 位 置 .目标 形状 
等 ,今后 将 根据 不 同 内 容 的 图 像 和 不 同 需求 的 分 割 任务 
探索 多 种 描述 算 子 的 表达 和 融合 ,以 实现 面向 应 用 需求 
的 更 好 分 割 。 
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